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ABSTRACT 

We present XMM-Newton observations and spectral fitting of two highly redshift, 
[OIII]-luminous, narrow-line radio galaxies, B3 0731+438 and 3C 257. Their X-ray 
continua are well fitted by a partial covering model with intrinsically unabsorbed 
and absorbed powerlaw components. The spectral models indicate that both objects 
harbour highly obscured nuclei, with Nh ~ 0.5 — 2 x lO^'^ cm^^. Correcting for this 



absorption we find large intrinsic luminosities in the range Lx 
Thus, both sources are type-2 quasars. 
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1 INTRODUCTION 

An important issue in AGN research is whether heavily ob- 
scured powerful AGNs (type-2 quasars) are common. The 
unification of AGNs (Antonucci 1993), UV spectral shapes 
(Maiolino et al. 2001) and the modeling of AGNs infrared 
spectral energy distributions (Pier & Krolik 1993) all sug- 
gest the presence of numerous type-2 quasars in the Uni- 
verse. These obscured AGN are thought to outnumber unob- 
scured objects by a factor of « 3:2 at redshifts out to z « 3.5 
(Willott et al. 2000; Fiore et al. 2001) and it is important 
to understand whether this ratio applies at higher redshifts 
and luminosities (Barcons et al. 1998; Almaini, Lawerence 
& Boyle 1999; Stern et al. 2002). Obscured AGN are also 
believed to be responsible for a significant fraction of the 
hard X-ray background (Comastri et al. 1995). Almost all 
of the soft X-ray background (< 2 keV) has been resolved 
into discrete sources, the majority of which are broad line 
(BL) type-1 QSO's. At higher energies, typical unobscured 
QSO spectra are too steep to explain the XRB and account 
for only ~ 20 % of the sources (Hasinger et al. 1998). Ob- 
scured AGNs provide a natural explanation, since photoelec- 
tric absorption allows only the hard X-rays to penetrate. 
This implies the majority of the energy density generated 
by accretion in the Universe takes place in obscured AGN. 
These results have been confirmed by data from two deep 
C/iandra observations (Brandt et al. 2001; Rosati et al. 2002) 
and recent analysis from an XMM-Newton observation of the 
Lockman Hole has also revealed similar results; Mainieri et 
al. (2002) found ~ 27 % of the Lockman Hole sample are 



Extremely Red Objects, the majority of which contain an 
obscured AGN. 

Recently a small number of type-2 quasars candidates 
have been identified, but these results are usually based 
on fimited data (Fabian et al. 2000; Willott et al. 2001; 
Fabian et al. 2001). The best published type-2 quasars, 
all obtained using the Chandra X-ray Observatory, are 
IRAS 09104+4109 at z = 0.44 (Iwasawa et al. 2001), orig- 
inally observed with BeppoSAX (Franceschini et al. 2000), 
CDF-S 202 at z = 3.7 (Norman et al. 2001), A18 at z = 1.467 
(Crawford et al. 2001), CX052, z = 3.288 (Stern et al. 2002) 
and B2 0902+343, z = 3.395 (Fabian, Crawford & Iwasawa 
2002). Correcting for absorption, these objects have X-ray 
luminosities of the order 10*^ erg s~^ and intrinsic absorp- 
tion Nh « 10^* cm-^ 

To identify likely type-2 quasars which are bright 
enough to provide reasonable X-ray spectra, we have ob- 
served two Narrow line radio galaxies (NLRGs) with XMM- 
Newton. The correlation between the mass of the central 
SMBHs and the galaxy spheroidal component (Magorrian et 
al. 1998) indicates that SMBHs are most hkely to reside in 
massive spheroids. Radio galaxies almost always host (form- 
ing) massive spheroids (van Bruegel et al. 1998), and NL- 
RGs lack broad components to their optical emission lines, 
so these galaxies are promising candidates for hosting ob- 
scured SMBHs (Barthel 1989). To identify such objects sev- 
eral parameters need to be measured; in particular the ab- 
sorbing column of the obscuring material and the intrin- 
sic X-ray luminosity. In the context of this paper Type-2 
quasars are defined to have an intrinsic, hard X-ray (2 - 
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Model 




r 


Nh 

(1022 cm-2) 


x2/d.0.f 


f-test 

% 


1 
2 
3 


B3 0731+438 

po X wabs 

po X wabs X zwabs 

wabs(C*po + zwabs x po) 


1 42+015 

1 QQ+0.54 
^•^'•^-0.10 

^•"1-0.37 


0.06 

Q o +2.5 
^•^ -2.6 


29.15/32 
26.11/31 
16.89/30 


93.30 
99.97 


1 
2 
3 


3C 257 

po x wabs 

po X wabs x zwabs 

wabs(C*po + zwabs x po) 


1 03+°" 
1 46+°-20 


0.04 
5.5 


7.15/12 
4.04/11 
3.96/11 


98.60 
98.60 



Table 1. The Models used to fit the 0.3 - 10 keV XMM-Newton data. For B3 0731+438 the scattering efficiency constant, C in Model 
3, was found to be 11 % ± 6 % and for 3C 257 it was fixed at 10 %. Also included in the table are the f-test results from comparing fits 
2 and 3 with fit 1. One sigma errors are quoted. 



10 keV) luminosity > 10** erg s^^ (we assume cosmological 
values Ho = 75 km Mpc~^ and qo = 0.5 throughout); 
Less luminous objects axe called Seyfert galaxies. The X-ray 
absorbing (gas and dust) column density of Seyfert 2s has 
Nh ^ 10^^ - 10^"^ cm"2 (Kigaliti, Maiohrio & Salvati 1999), 
and type-2 quasars are expected to have similar values. 

We present XMM-Newton observations and spectral fit- 
ting of two powerful FRII radio galaxies, B3 0731-^438 and 
3C 257. Both sources are high redshift NLRGs with pre- 
dicted X- ray luminosities well within the range of that ex- 
pected for type-2 quasars. 3C 257 is the highest redshift 
(z = 2.474) radio galaxy observed in the 3C catalog with 
a fairly smooth K band emission and morphology resem- 
bhng an elliptical galaxy (van Breugel 1998). B3 0731-1-438 
(z = 2.429) is a typical FRII radio galaxy with double-lobed, 
radio emitting hot spots and a central core (Caxilli et al. 
1997; Motohara et al. 2000). 

This paper is organized as follows; in Section 2 we dis- 
cuss sample selection and the use of optical emission line 
luminosities to determine the predicted X-ray luminosity of 
both sources. This is followed in Sections 3 and 4 by the 
XMM-Newton observations and data analysis. In Section 5 
we discuss the spectral fitting and in the final section we 
conclude our findings. 



nent (FWHMHa+[Nn] 800 km s~^) and 3C 257 has a some- 
what broader line (FWHMHa+[Nii] 1800 km s~^). Thus, for 
both objects their spectra indicate prominent emission from 
a narrow-line region while any nuclear emission is highly 
obscured. 

For B3 0731-1-438 the diagnostic emission-line ratios are 
consistent with the presence of a Seyfert-2 nucleus (Evans 
1998). A narrow-band optical image shows a large Lya emis- 
sion nebula (McCarthy 1991), suggestive of a hidden nucleus 
ionising the extended gas. This possibility is supported by 
a narrow-band 2.25 /xm image obtained using the Subaru 
telescope showing biconical clouds extending out to 37 kpc, 
radiating Ha-|-[NII] emission lines with a total Hq+[NII] 
flux of 3.5 x lO"^'^ erg s"^ cm"^ (Motohara et al. 2000), 
corresponding to a luminosity of 6.7 x 10*'^ erg s^^, (cor- 
recting to a cosmological scale of Ho = 75 km s~^ Mpc~^ 
and qo = 0.5). The Ha-|-[NII] ionization cone, aligned with 
the axis of the radio hot spots, implies the existence of a 
hidden, more luminous central source not observed in the 
optical, consistent with Keck spectropolarimetry data which 
suggest the nucleus is extinguished and we are only seeing a 
young stellar population in direct optical light (Vernet et al. 
2001). The combination of these data suggest B3 0731+438 
is a powerful, radio-loud AGN with a large intrinsic absorp- 
tion. Little is known about 3C 257. 



2 TARGET SELECTION 

There is an observed correlation between L[oiii] (A5007) 
and the X-ray luminosity Lx (2 - 10 keV) in radio 

loud AGNs (Mulchaey et al. 1994; Sambruna, Eracleous 
& Mushotzky 1999). This relation indicates NLRGs with 
L[oiii] > 5 x 10*'^ erg s^^ will have an X-ray luminosity 
Lx (2 - 10 keV) > 5 x 10** erg s"^ at the 90 % probabil- 
ity level. The relationship was used to select target NLRGs 
which are likely to host luminous AGNs and the two most lu- 
minous in [OIII] were chosen for observation. From the [OIII] 
luminosities (Evans 1998) the predicted X-ray luminosities 
for B3 07314-438 (L[oni] = 12.4 x 10*=^ erg s'^) and 3C 257 
(L[oiii] = 14-0 x 10*^ erg s"^) were estimated to be greater 
than 12 and 14 x 10** erg s^^ respectively at the 90 % prob- 
ability level. The width of the Ha line in the K-band spectra 
was determined for both sources (Evans 1998; Jackson and 
Rawlings 1997); B3 0731-1-438 shows no broad Ha compo- 



3 OBSERVATIONS AND DATA REDUCTION 

B3 0731-1-438 and 3C 257 were observed during revolutions 

343 (22nd October 2001) and 370 (15th December 2001) of 
XMM-Newton respectively. The observations of the EPIC 
MOS and PN detectors for both sources were in full window 
mode using a medium filter. The data were screened with 
the XMM Science Analysis Software (SAS) and pipeline pro- 
cessed. X-ray events corresponding to patterns 0-12 for the 
two MOS cameras and patterns 0-4, corresponding to single 
and double events for the PN were selected. 

We then extracted the spectra and lightcurves for both 
the source and background in each observation, separately 
for each EPIC detector. Circular source regions of 1 arcmin 
diameter for B3 0731-1-438 and 0.6 arcmin for 3C 257 were 
selected. The background regions were placed at positions 
close to the sources and any significant flaring in the source 
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and background regions were then removed using time se- 
lections in XMMSELECT. We then coadded the MOSl and 
M0S2 data into a single spectral file to maximize the signal- 
to-noise ratio and grouped the data into a minimum of 20 
counts per bin. 

The total effective observing time for B3 0731-1-438 is 
22.75 ksec for PN and 27.40 ksec for both MOSl and M0S2 
with a source count rate of 0.01 cts/s for PN and 0.0033 cts/s 
for the coadded MOS data. The effective observing time 
for 3C 257 is 25.73 ksec for PN and 30.15 ksec for both 
MOSl and MOS2, with a source count rate 0.0069 cts/s 
for PN and 0.0028 cts/s for the coadded MOS data. Finally, 
the background subtracted spectra were fitted using XSPEC 
version ll.l.Oae with the latest response matrices. 



4 SPECTRAL ANALYSIS 

The spectral analysis was carried out in the conventional 
way. The data were compared with a range of parametric 
models until the best fit (best reduced x^) obtained. 
Seperate PN and MOS fits were carried out as a standard 
check for the fit and it was found that the data produced 
similar results to that of the coadded PN and MOS data. 
Hereafter, we present fits only to the coadded data with one 
sigma errors quoted throughout. 

4.1 B3 0731+438 

Initially we fitted the coadded PN and MOS spectra 
(0.3 - 10 keV) with a single powerlaw and neutral ab- 
sorption (Galactic line-of-sight absorption) model; Model 

I in Table 1. Table 1 also includes values for the pow- 
erlaw photon index (F), the associated intrinsic column 
density (Nh), where appropriate, and values for ''■nd 
number of degrees of freedom in the fit. Using a galac- 
tic column density Nh = 6.2 x 10^° cm"^, Model 1 
gives F = 1.4 ± 0.15 {x^/d.o.i = 29.15/32), suggesting 
some intrinsic absorption. The unabsorbed fiux and lu- 
minosity (2 - 10 keV) were found from the fit to be 

3.2 X 10"" erg s"^ cm~^ and 2.8 x 10** erg s"^ respec- 
tively. 

The fit is improved (93 % significance) when intrin- 
sic absorption is included in the spectral model; Model 2. 
This produces a steeper underlying powerlaw, F ~ 1.8 and 
Nh « 3.3 X 10^^ cm"^ (x^/d.o.i = 26.11/31). This resuh 
suggests the source has a possible intrinisically scattered 
component, typical of a type-2 AGN. Therefore a partial 
covering model was fitted; Model 3, consisting of intrinsically 
absorbed and unabsorbed (scattered) powerlaw components. 
Model 3 provides an improved spectral fit for B3 073H-438 
at greater than 99.9 % significance, with F = 2.6 ± 0.4 and 
Nh = 1.9 X 10^^ cm-2 ± 0.7 {x^/d.o.i = 16.89/30). The un- 
absorbed X-ray luminosity (2 - 10 keV) is 9.1 x 10** erg s~^. 
Figure 1 displays the best fit (Model 3) for B3 0731-1-438; the 
partial covering model. Figure 2 illustrates a contour plot of 
Model 3, showing the column density and X-ray luminos- 
ity distribution. The scattering efficiency was found to be 

II % ± 6 %. Constraining the scattering efficiency to 10 % 
we find the results are consistent with the unconstrained fit. 
These results are consistent with the low scattering efficien- 
cies for Seyfert 2s (Lumsden et al. 2001). 
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Figure 1. The 0.3 - 10 keV EPIC coadded MOS (MOSl -|- 
MOS2) and PN spectrum of B3 0731+438 fitted with a partial 
covering model (Model 3 in Table 1). The lower plot shows the 
ratio of the data to the fitted model. For clarity, the data have 
been rebinned to 30 counts per bin. 
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Figure 2. A contour plot of column density against X-ray lu- 
minosity for B 0731+438 from Model 3. The confidence intervals 
shown are for 68, 90 and 99 %. 



The spectra of B3 0731+438 may be consistent with a 
refiection dominated component, but fitting the data with 
the refiection model PEXRAV (Magdziarz & Zdziarski 1995) 
requires a very steep photon index of F « 3.3 ± 0.2 and 
no direct view of the incident continuum. If we assume the 
refiected component is « 5 %, the implied upper limit on the 
intrinsic X-ray luminosity (2 - 10 keV), is 5.6 x 10*^ erg s~^ 
if the planar refiection subtends 2-k sr at the illuminating 
source. Fitting the data with a broad (1 keV) FeKa emission 
line, we find a 90 % confidence upper limit on the intensity 
to be 2.7 X 10~® photon s"'^ cm~^, corresponding to a rest- 
frame equivalent width of 263 eV. This equivalent width is 
low, particularly for a type-2 object (Nandra et al. 1997; 
Turner et al. 1997). 



4.2 3C 257 

The analysis of the PN data for 3C 257 was complicated by 
a residual spectral feature at « 1.5 keV due to incomplete 
subtraction of an internal detector background aluminium 
K-alpha line (Lumb et al. 2002). This line was removed from 
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channel energy (keV) 

Figure 3. The 0.3 - 10 keV EPIC coadded MOS and PN spec- 
trum of 3C 257 fitted with a partial covering model (Model 3, 
Table 1) binned to 20 counts per bin. The fit excludes the region 
from 1.4 — 1.6 keV (see text). 



the spectral fitting by ignoring energies in the range 1.4 - 
1.6 keV for the 3C 257 PN data during the analysis. 

The models used for B3 0731+438 were again used for 
3C 257. Model 1, the single powerlaw and neutral absorp- 
tion model, using a Galactic value of Nh = 4.4 x 10^" cm~^, 
produced a very hard photon index F = 1.0 ± 0.1 
(x^/d.o.f = 7.15/12). The flux and unabsorbed lu- 
minosity (2 - 10 keV) obtained from the fit were 
2.6 X 10~^* erg s~^ cm~^ and 1.5 x 10'*"' erg s~^ respec- 
tively. As in the case of B3 07314-438 the fit for 3C 257 is 
improved, with a 90 % significance, when intrinsic absorp- 
tion is included in the spectral model; Model 2. This again 
produced a steeper underlying powerlaw with F ~ 1.5 and 
Nh 5.0 X 10^^ cm"^ ix^/d.o.i = 4.04/11). 

We then used the partial covering model, Model 3, 
but it is not possible to distinguish between this model 
and Model 2. Since it was difficult to constrain the 
fit, we fixed the scattering efficiency to 10 %, consis- 
tent with that of B3 0731-1-438. We found F « 1.4 and 
Nh ~ 5.5 X lO^'^ cm~^ with a corrected, unabsorbed X- 
ray luminosity (2 - 10 keV) = 2.6 x lO*"* erg s~*. Figure 3 
displays the spectrum of source 3C 257, fitted using Model 
3. Both models 2 and 3 suggest there is significant intrinsic 
absorption in 3C 257. 

Fitting the data with the refiection model PEXRAV 
requires a steep photon index F « 2.7 with no direct view 
of the incident continuum. Making the same assumptions as 
for B3 0731+438 , the implied upper limit on the intrinsic 



X-ray luminosity (2 - 10 keV), is 3.0 x 10*'^ erg s"^ 



5 DISCUSSION 

The best-fit partial-scattering model for B3 0731+438 gives 
an intrinsic absorbing column of Nh ~ 1.9 x 10'^'' cm~'^ 
which is comparable to the column density suggested by 
Risaliti, Maiolino & Salvati, (1999) for type-2 Seyferts of 
Nh > 10'^'^ cm~^. This model also gave a photon index of 
F PS 2.6. This value is steep, but is comparable to the typical 
quasar value of F = 1.9, at 90 % significance. Some luminous 
quasars have steep photon indices (e.g. F ~ 2 for the quasar 
PDS 456 (Reeves et al. 1999)). For 3C 257 it was difiicuh to 



constrain the fit as a result of the limited number of photons. 
However, the spectral fits are consistent with this source also 
having a large X-ray luminosity and intrinsic absorption. 

Both sources appear less luminous than predicted based 
on their [OUT] optical fiuxes, but are still within the quasar 
range. The optical to X-ray spectral index Uox was deter- 
mined for B3 0731+438 using IR spectroscopy from Bales 
& Rawlings (1993) and Evans (1998). Assuming an optical 
extinction factor of Av = 1.9 (Motohara et al. 2000) and 
using an absorption corrected X-ray fiux, we calculate 
Qox = - 1.87 ± 0.2, consistent with sources studied by Bech- 
told et al. (2002). 

For B3 0731+438 we can compare the ionizing photon 
number derived by Motohara et al. (2000) with the observed 
X-ray continuum fiux. Integrating a powerlaw of photon in- 
dex 2.5, consistent with that used by Motohara et al. (2000), 
between 13.67 eV and infinity, normalised to the unabsorbed 
X-ray continuum, we derive a value of 1.7 x 10^^ pho- 
tons s~*. Repeating this calculation with our derived value 
of Qox, extrapolating to 13.67 eV using the SED of Zheng et 
al. (1997), we derive a value of 1.0 x 10^^ photons s~^ . These 
photon numbers are consistent with the result obtained by 
Motohara et al. (2000), of > 2.1 x 10^** photons s"* based 
on the Ha luminosity, assuming case B and a unity covering 
factor in the Ha cone (corrected to a cosmological scale of 
Ho = 75 km s~* Mpc~* and qo = 0.5). We can also com- 
pare the X-ray and Ha luminosities assuming the relation 
that the intrinsic X-ray luminosity (2 - 10 keV) is approx- 
imately ten times that of the Ha luminosity (Ward et al. 
1988; Ho et al. 2001). Again correcting to a cosmological 
scale of Ho = 75 km s"'^ Mpc"'^ and qo ~ 0.5, Motohara 
et al. (2000) found a Ha luminosity for B3 0731+438 of 
6.7 X 10'*'' erg s^*, giving a predicted X-ray luminosity of 
6.7 X lO'** erg s^*. This is in good agreement with our value 
obtained for B3 0731+438 of 9.1 x 10*** erg s-\ These com- 
parisons suggest the X-ray emission is fairly isotropic. 



6 CONCLUSIONS 

In this paper we have presented XMM-Newton spectral anal- 
ysis for the two NLRGs B3 0731+438 and 3C 257. The anal- 
ysis of B3 0731+438 from fitting the spectra with a partial 
covering model indicates a spectral index of F = 2.6 and a 
large column density Nh « 2 x 10^^ cm~^. As a result of 
few counts, it was difficult to constrain the fit for 3C 257. 
However, the results suggest a large column density with a 
hard photon index, in agreement with an obscured source. 
Both objects have large absorption corrected X-ray lumi- 
nosities (2 - 10 keV of 9 X 10** erg s"* for B3 0731+438 
and 2.6 x 10** erg s"* for 3C 257. 

In summary, the analysis has shown that both objects 
have intrinsic column densities and large X-ray luminosities 
well within the type-2 quasar range. We therefore conclude 
B3 0731+438 and 3C 257 are examples of type-2 quasar and 
currently provide the best X-ray spectra obtained of such 
objects with XMM-Newton. 
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